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Abstract 

A previous proof of non existence of tokamak equilibria with purely 
poloidal flow within macroscopic theory [Throumoulopoulos, Weitzner, 
Tasso, Physics of Plasmas 13, 122501 (2006)] motivated this micro- 
scopic analysis near magnetic axis for toroidal and " straight" tokamak 
plasmas. Despite the new exact solutions of Vlasov's equation found 
here, the structure of macroscopic flows remains elusive. 



1 



1 Introduction 



Some time ago (see [H[2]), it was possible to prove non existence of tokamak 
equilibria with purely poloidal incompressible flow. Recently, an extension to 
compressible plasmas appeared in Ref.[3J including Hall term and pressure 
anisotropy. The proof for the incompressible case given in Refs.jU [2] was 
global while the recent proof [3] is limited to the neighbouring of the magnetic 
axis through a kind of Mercier expansion. 

This last result motivated the idea to extend the analysis to Vlasov- 
Maxwell equations examined near axis. An important ingredient is to write 
the Vlasov equation in cylindrical coordinates in a tokamak geometry, which 
simplifies the subsequent analysis. We use for that purpose the calculation 
done in an old ICTP report [I] where the Vlasov equation is written in 
arbitrary orthogonal coordinates. 

In Section 2 the expression of the Vlasov equation is obtained in toroidal 
geometry. In Section 3 the ODEs of the characteristics are derived while 
Section 4 is devoted to "straight tokamaks" and Section 5 to discussion and 
conclusions. 

2 Vlasov equation in orthogonal coordinates 

As explained in Ref.Jl] we consider a general system of orthogonal coordi- 
nates x 1 , x 2 , x 3 with the metric ds 2 = gn^dx 1 ) 2 + g22(dx 2 ) 2 + g^dx 3 ) 2 and 
unit vectors = where i goes from 1 to 3. The velocity vector of a 
"microscopic" fluid element is then projected on the unit vectors e.; as 



where the components v l are independent upon space variables. The total 
derivative of v is 



where E and B are the electric and magnetic fields consistent with Maxwell 
equations and the charge to mass ratio — is set to one. Projecting Eq.(2) on 
the unit vectors we obtain 



v = v'e. 



(1) 



— + v-Vv = E + vxB, 



(2) 



dv l 
~dt 



e r (E + vxB) + e r vxVxv. 



(3) 
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Finally, the Vlasov equation in orthogonal coordinates is given by 



| + v.V/ + e r (E + vxB)| + (e ! .vxVxv)| = 0, (4) 

where / is a function of the x l , v l and time while v is given by Eq.(l). For 
more details see Ref. [4]. / stays here for the ion distribution while the 
distribution function for the electrons is governed by an equation similar to 
Eq.(4). 

Let us now specialize on cylindrical coordinates x 1 = r, x 2 = (ft, x 3 = z. 
Then V x e< = for % = 1 and 3 and V x e 2 = ei x V</>. If we replace the 
indices 1, 2, 3 by r, (ft, z we have V x v = v^e r x V0 and 

v r v% <y) 2 e r 
v x V x v = . (5) 

So the last term of Eq.(4) becomes [^r~§^ ~^~^"§^\- Setting B = near 
axis and |£ = for steady state, Eq.(4) reads 



^ df [v z I - (v ^) 2 ] df | v r I df v r v^> df _ Q 

Assuming V/ = V$ = on axis the final equation to solve is 

- [v z I ~ (^) 2 # - + v r I^- = 0. (7) 

L V I iQ v r Q v 4, g v z \ I 

3 ODEs for characteristics 

Let us start with the simpler case 1 = 0, then the characteristics of Eq.(7) 
are given by the solution of 

dv r dv* 

(8) 



f v 4>)2 v r y 



whose solution is (t/) 2 + «) 2 = C. Since / = f(C,v z ) = f[((v r ) 2 + (v^) 2 ),v z ] 
on axis we obtain for the toroidal flow 



J v^fd 3 ^ = 0, (9) 



which means zero toroidal flow on axis. 

For 1^0 the charasteristics are given by 

dv r dv^ dv z 

v zj _ ( v <t>y ~ v r v 4> ~ v rj- y > 

The last equality delivers C\ = v z + I In \ v^\, the second characteristic being 
the particle energy C2 = {v r ) 2 + (v^) 2 + (v z ) 2 . C\ is "antisymmetric" in v z 
but symmetric in v^, which leads to 

/ ^/(Ci, C 2 )d\ = 0, J v z f(C^ C 2 )d 3 v ^ 0. (11) 

It means that the ^-flow is zero while the unphysical z-flow is finite. This is 
obviously not acceptable. 



4 "Straight" tokamaks 

The straight tokamaks do have magnetohydrodynamic solutions with purely 
poloidal flow as known from previous work [5]. For the purpose of a micro- 
scopic theory the appropriate coordinate system is the cartesian one x 1 = x, 
x 2 = y, x 3 = z so that the toroidal angular coordinate is replaced by y and 
the toroidal field I by B y . Since Vxe^ vanishes for all i, the term v x V x v 
in Eq.(4) disappears. 

For the steady state with finite B y , Eq.(7) is replaced by 

-^ + ^ = °- (12) 



whose characteristic is given by 

dv x dv z 



V' 



(13) 



The solution of Eq.(13) is C = (v x ) 2 + (v z ) 2 , which leads to / = f((v x ) 2 + 
(v z ) 2 , v y ). Purely poloidal flows are possible, which is consistent with Ref. [5]. 
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5 Discussion and Conclusions 



The result of section 3 obliges us to change the assumptions leading from 
Eq.(6) to Eq.(7) i.e. V/ 7^ instead of zero on the magnetic axis. The 
special canonical 0-momentum solution is of that kind, and leads naturally 
to toroidal flows but no poloidal flows. However, a comprehensive discussion 
of the problem cannot be done since the complete set of characteristics of 
Eq.(6) is not known. 

Finally, though we know from section 3 that / must be a function of C\ 
and C2, we could, in addition, choose / to have different values for different 
signs of, for instance, rr. A known example of that kind of solutions is the 
case of BGK waves [5J, in which the "free particles" have different distribu- 
tions for different signs of their velocities. See also Ref. [7J for a quasi-neutral 
treatment. Though toroidal flows can then be constructed, physical con- 
straints like isotropy of the pressure tensor or constraints on other moments 
or geometrical symmetries and, ultimately, collisions could exclude such so- 
lutions. Again we are led to look for the general solution of Eq.(6) with 
V/ 7^ on axis in order to discuss the structure of the macroscopic flows. 
Unfortunately, as mentioned before, the answer to this problem is quite un- 
certain. 
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